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The dynamics of an elongated attractive Bose-Einstein condensate in an axisymmetric harmonic 
trap is studied. It is shown that density fringes caused by self-interference of the condensate order 
parameter seed modulational instability. The latter has novel features in contradistinction to the 
usual homogeneous case known from nonlinear fiber optics. Several open questions in the interpre- 
tation of the recent creation of the first matter-wave bright soliton train [Strecker et al. Nature 
417 150 (2002)] are addressed. It is shown that primary transverse collapse, followed by secondary 
collapse induced by soliton-soliton interactions, produce bursts of hot atoms at different time scales. 
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Solitons are ubiquitous in nature, appearing in systems 
as diverse as DNA, shallow water, and laser light pulses in 
nonlinear fiber optics 0. Modulational instability (MI), 
an effect well known in the latter context, is the process 
by which a constant-wave background becomes unstable 
to sinusoidal modulations due to the presence of a focus- 
ing nonlinearity, leading to a pulsed wave, called a bright 
soliton train [2|. Recently, the first matter- wave bright 
solitons were created from attractive Bose-Einstein con- 
densates (BEC's) in elongated harmonic traps 0,0. In 
the latter experiment a train of from four to ten soli- 
tons was formed; however, unlike the solitons resulting 
from MI of a uniform initial wavefunction in a constant 
potential in one dimension , these solitons were appar- 
ently stabilized by repulsive soliton-soliton interactions. 
Moreover, of the ~ 3 x 10 5 initial atoms present when the 
focusing nonlinearity was turned on, only ~ 10% survived 
to form the soliton final state, implying massive collapse. 
It was suggested 00 that MI, in this new context, cre- 
ates such a repulsively interacting soliton train. 

In the following, the MI of a non-uniform initial state 
in the presence of a harmonic potential is studied both 
analytically and numerically in the context of the mean 
field of the BEC. The conditions by which the transverse 
dimensions may bring about primary collapse, either of 
the entire condensate in the radial direction, or of indi- 
vidual solitons as they are formed by MI, as well as sec- 
ondary collapse induced by soliton-soliton interactions, 
are detailed. As a matter-wave bright soliton train is 
in fact a train of self-contained BEC's easily guided and 
manipulated by electromagnetic fields, it is important to 
understand how it is made and how it may be kept stable. 
Moreover, bright solitons, which have already revolution- 
ized the communications industry 0, are, in the context 
of the BEC, predicted to have applications in atom in- 



terferometry and quantum frequency standards jaj- 

The results of our analysis differ strongly from those 
of previous studies HI ■ Firstly, it is shown that self- 
interference of the order parameter seeds MI at different 
times and with a steadily decreasing wavelength. This is 
true for any non-uniform initial density with the excep- 
tion of a Gaussian. Secondly, it is shown that the ensuing 
solitons may have any relative phase, and may there- 
fore interact both repulsively and attractively. Thirdly, 
it is shown that both primary and secondary collapse oc- 
curs. This leads to the testable experimental prediction 
of bursts of hot atoms emitted at different time scales and 
to the conjecture that trajectories of remaining solitons 
are stabilized by selection, as detailed below. 

The 3D Nonlinear Schrodingcr equation (NLS) or 
Gross-Pitacvskii equation which describes the mean field 
of the BEC is written as 

[-fi 2 V 2 /2m + .g TV |*| 2 + V(r)]$ = ihd t V , (1) 

where V(r) = m(UpP 2 + uj 2 z 2 )/2, g = 4nK a s /m, a s is 
the s-wave scattering length, m is the atomic mass, TV 
is the number of condensed atoms, the condensate or- 
der parameter "3/ = 4'(r, t) has been normalized to one, 
and axisymmetric confinement has been assumed. Note 
that for negative scattering length, or attractive nonlin- 
earity, solutions are liable to collapse in certain param- 
eter regimes, as shall be discussed below. In the case of 
strongly anisotropic confinement, one may adiabatically 
separate the slow longitudinal from the fast transverse 
degrees of freedom. The adiabatically varying transverse 
state obeys a 2D NLS which shows an instability towards 
collapse. The criterion for stability found by numerical 
integration of the 2D NLS is 
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where NUMz, t)\ 2 is the local axial line density of the con- 
densate |l5j. If adiabaticity is violated, collapse can also 
happen at weaker nonlinearity due to transverse oscilla- 
tions on a time scale tt/uj p |8|. When the longitudinal 
dynamics is significantly slower than this time scale, the 
adiabatic separation is valid. If, additionally, the trans- 
verse nonlinearity is weak, i.e., \n\ <C 1, the longitudinal 
equation reduces to the quasi-lD NLS 

[-h 2 d 2 /2m- 



sin(fc_z + 8 - § ) 
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where (?id = 2a s uj p T% is the quasi-lD coupling con- 
stant 0, provided l p S> \a\ [T^j, with l p = ^/hj (moj p ). 

In order to understand the mechanism of MI for a non- 
uniform initial density profile and in the presence of a 
non-constant potential, it is necessary to briefly review 
MI in the uniform case, which is well known from fiber 
optics 0]. A linear response analysis reveals that, for 
attractive nonlinearity, a small sinusoidal modulation of a 
uniform state tpo with wavenumber k grows exponentially 
at a rate 7 given by 

2 2AH 5lD ||Vo| 2 l 2 , l^'ol 4 ^ 21 
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The maximum growth rate 7, ng = 
obtained at wavenumber k. 
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\a s \N\^ \ 2 is 
where £ = 



l p / y/4\a s \ N \ipo\ 2 is the effective ID healing length of 
the condensate [3~1] . Note that iV|-0o| 2 is the line den- 
sity of the condensate. Growth occurs only if j 2 > 0, 
which implies < k < fc max = y/2k mg . This means that 
nonlinear focusing can only be seeded by modulations of 
sufficiently long wavelength and is fastest at the length 
scale of 27r£. 

For non-uniform initial density profiles, self- 
interference fringes in the order parameter can provide 
the necessary seed. To understand this phenomenon, 
it is instructive to first consider the development of 
the wavefunction for the ideal gas (a s = 0). As a 
model, we shall consider the case of an initial state 
ip{z, 0) which is a rectangular function of height ipo 
and width L centered at the origin. Recall that the 
Feynman propagator which determines the evolution of 
the wavefunction in the linear Schrodingcr equation is 
defined by ip(z,t) = J dz'G(z,t; z\0)ip(z',0). In the 
case of a harmonic potential, G may be determined 
exactly by semi-classical methods to be 01 03 , 

^ cxp {i(z 2 - 2zz'/cost + z' 2 )/(2^tanr)} 
l z ^/2iti \ sin"r| 



where r = u> z t and l z = y/h/(mu> z ) is the axial harmonic 
oscillator length. This gives an analytic expression for 
the evolution of the order parameter, which may be more 
easily understood in the limit ui z t <C 1: 
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where 9 = Arg(i/)). Figures ^b) andJS^b) show the den- 
sity and phase structure given by Eqs. (101811 . respectively. 
Note that the above expansions converge very slowly as 
\z\ — > L/2 where the initial wavefunction is dicontinuous. 

The above expressions for the evolution of the phase 
and density of the wavefunction may be understood as 
follows. All trigonometric terms represent quantum self- 
interference and produce fringes which seed MI. The 
wavenumber k± is dependent on both time and position; 
the wavelength is longer and the amplitude of oscillations 
higher at the edges of the box than in the center, and the 
overall wavelength increases as a function of time. As 
these interference terms are to linear order independent 
of lo z (since l z ui z = h/m), they must result solely from 
the non-uniform initial density profile. One may also un- 
derstand this result by observing that Eq. JSJ is, to order 
u> z t, identical with the free space Feynman propagator in 
one dimension. In contrast, the non-trigonometric terms 
are caused by the harmonic potential: in Eq. © they 
lead to a monotonic increase in the mean density, where 
the order (w 2 i) 2 term is independent of position; and in 
Eq. J7J) they develop the overall phase profile harmoni- 
cally. The development of the density and phase of the 
wavefunction for the geometry of Ref. QJ of L ~ 10 l z is 
shown in Figs. 2] and 

MI of a non-uniform initial state, both with and with- 
out an additional harmonic potential, may be understood 
in light of the above considerations. In Fig. [3] is shown 
the numerical evolution of the axial line density N\?jj\ 2 in 
Eq. © for the parameters of Strcckcr et al. |j] and our 
model of an initial square function. The fringes produced 
by self-interference of the condensate order parameter tp 
attain the critical wavelength for MI of y/2n£ first at 
the edges of the cloud, as shown in Figs. and In 
Fig. |21 one notes that the solitons are therefore formed 
last in the center of the trap. In this latter stage they 
are also formed closer together, as, due both to the har- 
monic potential (see Eq. (7J) and the global focusing ef- 
fect of the nonlinearity, the axial density has increased, 



3 



0.4 
h 
0.3 

0.2 
0.1 



(a) 



_2D_cpllap_se Limit _ 



0.1 



2.5 


z/ l z 5 




A/\/\J 


(b) 





2.5 



z/l z 5 



FIG. 1: Shown is the right half of the axial density of a 
BEC in a harmonic potential, as determined by (a) numerical 
solution of the nonlinear Schrodinger equation, and (b) an- 
alytical solution of the linear Schrodinger equation (dashed: 
exact, solid: approximate), at t = 0.016 x 2ir/u> p , after ID 
evolution from an initial box state of length L/l z = 10 cen- 
tered at the origin. The formation of solitons first at the edge 
of the condensate in (a) results from the longer wavelength 
of the linear fringes, as seen in (b), since MI takes place at 
a wavelength ~ 2-7r£ = 0A3l z in this simulation, using the 
parameters of Strecker et al. Fil l 1(1 . 



and £ oc IV'I -1 - Note that although the early phase struc- 
ture in Fig |21 depends only weakly on the harmonic po- 
tential, in later stages of evolution the trap affects the 
phase strongly. 

One may ask if the seed fringes are caused by the sharp- 
ness of the square function edges in our model. In the 
particular case uj z t = 7r/4, Eq. JSJ takes the form of a 
Fourier transform, tp(z 7 t) oc J dz' exp(— i zz' /l^)ip(z' , 0). 
Thus it is immediately apparent that all initial wave- 
functions with the exception of a Gaussian must develop 
fringes. In order to test the full time development for an 
experimentally relevant initial density, a Thomas-Fermi 
profile of harmonic trapping with ip oc yj\ — (2z/L) 2 
for \z\ < L/2 and i/j = for \z\ > L/2 was also studied 
and was found to have qualitatively the same evolution 
in the linear Schrodinger equation, i.e., fringes of longer 
wavelength towards the outer edges of the density pro- 
file and an overall increase in wavelength with time. The 
nonlinear evolution also shows the same qualitative be- 
havior as Fig. |3 with solitons forming spontaneously first 
near the edges and later toward the center. 

Besides the initial linear density fringes and nonlinear 
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FIG. 2: Shown is the phase of the order parameter corre- 
sponding to Fig. H^a) and (b). The phase difference A<f> be- 
tween solitons may be seen by comparing to the large density 
peaks of Fig^a): A</> is highly sensitive to the nonlinear- 
ity, and may take on arbitrary values between and 2n, in 
contrast to the usual case of MI in fiber optics, where it is 
restricted to 0. At later stages in the simulation the relative 
phases drift as the solitons decouple. 



MI, the condensate can undergo primary transverse col- 
lapse on a time scale tt/lu p , according to Eq. (J2J). In the 
Strecker experiment, a rough estimate assuming a con- 
stant initial density yields r\ 6.3. In fact, r\ will exceed 
this value locally due to the slanted initial profile pro- 
duced by creating the BEC on one side of the trap £| 
and transverse collapse is likely to occur on the higher 
density side of the profile before MI develops. However, 
even for a uniform initial state, the critical value of r\ is 
reached during soliton formation as indicated in Fig. J5J. 
Note that, the center of mass of the BEC, created at the 
side of the trap, trivially oscillates with the harmonic 
oscillator frequency and decouples completely [l3T | from 
the relative motion, the latter of which is analyzed in this 
work. 

After the soliton train has spontaneously formed, even 
if transverse collapse is avoided, individual solitons can 
undergo primary 3D collapse, which dominates over 
transverse instabilities. The static condition for three- 
dimensional soliton stability in the absence of longitudi- 
nal confinement is given numerically by Q 



N ciit \a s \/l p = 0.627. 



(9) 



Dynamical effects, as for example breathing of the soli- 
ton, can lower this condition further. In light of a recent 
experiment on a single matter- wave bright soliton 
it is anticipated that such a three dimensional collapse 
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FIG. 3: Shown is the ID evolution of the axial density of 
an elongated Bose- Einstein condensate when the interactions 
are suddenly tuned negative, starting from a box-like initial 
state. Here the parameters of the Strecker experiment Q are 
used. At early times, small fringes develop as the time-slice of 
Fig. clearly shows; once the MI is seeded, the nonlinearity 
focuses the fringes to produce solitons, the leading edge of 
which is seen at the bottom of the figure as a wide triangle 
shape. Lines visible above this edge are solitons. Their com- 
plicated interaction pattern clearly shows many examples of 
both attractive and repulsive soliton interactions. 



will eject the excess number of atoms, leaving behind a 
solitonic core which is near the limit of Eq. ©. In the 
simulation of Figs. 11131 the 26 solitons formed carry an 
average of « 1.9iV cr i t atoms where 7V C rit « 6100 for the 
parameters of the Strecker experiment. These results are 
supported by a coarse estimate based on Eq. (J1J, which 
predicts formation of 22 solitons, and by simulations of 
the 3D NLS analogous to the one of Figs. 11131 which 
shows primary collapse immediately after formation of 
the first solitons (not shown). 

Finally, the coherent overlap of solitons, not themselves 
subject to two or three dimensional collapse, can cause 
secondary collapse. The dynamics of binary soliton in- 
teractions are determined by their relative phase A<p and 
amplitude A A Q. For -7r/2 < A<j> < it/ 2 they attract 
each other and exchange mass. Consequently, they can 
undergo collapse if, during their interaction, they tem- 
porarily violate Eq. ©. For ir/2 < A(f> < 3ir/2 and 
AA = the interaction is repulsive and there is no mass 
exchange; the solitons are stable against secondary col- 
lapse induced by interactions. However, for A A ^ the 
relative phase cycles periodically which, if it crosses over 
into the regime — ir/2 < A(j> < ir/2, can again cause 
collapse. Based on the above and the experiment of 
Ref. [3j, we conjecture that, in the Strecker experiment, 
collisions of solitons close to iV cr it led, by several stages 
of secondary collapse, to a final configuration of a small 
number of solitons with stable trajectories. We empha- 
size that these results differ from those found in Ref. |5j, 



where it was suggested that the phase difference A<f> be- 
tween adjacent solitons originates from quantum fluctua- 
tions and is restricted to values close to ir, thus stabilizing 
their trajectories. We have shown that MI of a nonuni- 
form initial state may produce arbitrary A<j) already in a 
mean-field model. We conjecture that it is the secondary 
collapse processes which determine a stable trajectory of 
solitons in the final state with tt/2 < A</> < 37r/2. 

The 2D and 3D collapse processes of Eqs. and (jSJ 
occur on different time scales. It is likely that both pro- 
cesses contributed to the ~ 90% initial atom loss in the 
Strecker experiment; a measurement of the time depen- 
dence of the atomic mass ejected from the trap during 
the initial stages of soliton train formation may be able 
to determine if one or both were dominant. For the ex- 
perimental parameters, based on the collapse conditions 
and Fig. primary transverse collapse occurs on a time 
scale of < 1.7 ms, primary 3D collapse occurs between 
1.5 and 10 ms, and secondary collapse occurs at > 10 ms. 

In conclusion, the phenomenon of MI has been ana- 
lyzed in the case of a non-uniform initial state in the 
nonlinear Schrodingcr equation with a harmonic poten- 
tial. It was shown that linear density fringes, the length 
scale of which depends both on space and time, seed the 
nonlinear instability. Primary transverse collapse, pri- 
mary three-dimensional soliton collapse, and secondary 
three-dimensional collapse due to soliton binary inter- 
actions were discussed, and all were predicted to have 
played a role in the Strecker experiment 0. 
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